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Abstract 

Using a quark model, we study the effect of charge symmetry breaking on the 
valence quark distributions of the nucleon. The effect due to quark mass differences 
and the Coulomb interaction of the electrically charged quarks is calculated and, in 
contrast to recent claims, found to be small. In addition, we investigate the effect 
of charge symmetry breaking in the confining interaction, and in the perturbative 
evolution equations used to relate the quark model distributions to experiment. We 
find that both these effects are small, and that the strong charge symmetry breaking 
effect included in the scalar confining interactions may be distinguishable from that 
generated by quark mass differences. 







1 Introduction 



The study of charge symmetry and its breaking is almost as old as nuclear physics itself. From 
the earliest days of isospin0, to modern attempts to understand small effects in the nucleon- 
nucleon interaction |fj], the study of this symmetry has provided a rare window into the 
non-perturbative dynamics of low energy hadronic phenomena. The interactions responsible 
for charge symmetry breaking (CSB) are largely understood and relatively weak, so that 
the study of CSB provides a sensitive filter with which to test the hadronic wavefunctions 
of nuclei and nucleons. In the case of QCD, where the interactions that bind quarks into 
hadrons are only understood schematically and the theoretical landscape is cluttered with 
different phenomenological models, such probes may prove especially valuable. 

In this paper, we study the effect of the breaking of charge symmetry on the valence 
quark distributions of the nucleon. Our interest in this topic is generated by the observation 
of unexpected effects in both the sea|| and spin-dependent 0] quark distributions of the 
nucleon, and by the possibility that large CSB effects in the valence distributions of the 
nucleon may play a non-negligible role in the extraction of sin 2 9yy from v — N data||. 
Recent calculations have claimed that charge symmetry breaking may be greatly enhanced 
by kinematic effects associated with diquarks in the nucleon wavef unction^], and that the 
resulting CSB may be large enough to be observed directly. Our calculations provide an 
important check on the model dependence of these results, as well as investigating several 
sources of CSB not considered in Ref. 6. 

The paper is structured as follows: In Section 2, the method we use to extract quark 
distributions from quark models is reviewed, with particular emphasis on those details of the 
calculation that will be affected by charge symmetry breaking. In the next section, the shift 
in the valence quark distributions due to the breaking of charge symmetry by quark masses, 
and by Coulomb effects is calculated for both the minority and majority quark distributions 
in the nucleon. The effect of the difference between the neutron and proton mass is also 
discussed and, after concluding that it contains no physics, discarded. The fourth section 
is devoted to a discussion of the possibilty of charge symmetry breaking in the confining 
interaction itself, and the possibility that such effects may be distinguishable from ordinary 
quark mass effects on the valence distributions, while in the fifth we calculate the charge 
symmetry breaking effect due to the perturbative evolution of the valence distributions that 
is required if the quark model is to make contact with high energy data. The final section 
compares our results to those of references 5 and 6, and discusses the prospects for measuring 
CSB in the valence distributions directly. 

2 Valence Quark Distributions in the Los Alamos Model 
Potential (LAMP). 

Quark models were originally constructed to provide a description of low energy hadronic 
data using only effective interactions between valence quarks. In order to make sensible 
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calculations of valence quark distributions, the stark simplicity of quark models must be 
reconciled with the richer vision afforded by data from higher energies, where nucleons are 
composed not only of valence quarks, but also of sea quarks and gluons. Remarkably, 
these two very different pictures may be accomodated by QCD via the renormalization 
group [0]. Quark models may be interpreted as representations of QCD at an intermediate 
renormalization scale, /Uq M , large enough so that the quark substructure of the nucleon is 
revealed, but small enough that the sea quarks and glue are almost entirely absorbed by a 
redefinition of the valence quarks. Parton distributions at this intermediate scale may be 
calculated in terms of the quark model, and then evolved to higher energies using perturbative 
QCD[[8] and compared to data. 

The first step in this procedure is to calculate the quark distributions at the quark model 
scale. For unpolarized scattering, the relevant matrix elements are given by|| 
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where q + = —Mx/y/2 (with x = Xsj the Bjorken scaling variable), ^(^i) are ne ld operators 
for quarks of flavor i, 7 + is a Dirac gamma matrix, and the subscript LC indicates a light 
cone condition on £, namely that £ + = £j_ = 0. 

For a given model, valence quark distributions can be obtained via a number of different 

The approach we adopt consists of 
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prescriptions that have arisen in the literature | 
a straightforward evaluation of the matrix elements of Eq. 1 in a Peierls-Yoccoz projected 
momentum eigenstate, assuming that the time dependence of the field operator is dominated 
by the lowest eigenvalue of the Dirac equation used to obtain the wavefunctions of the struck 
quark. The details of this procedure are described in Ref. 9, where the valence quark 
distributions are shown to be given by 



xq v {x) 
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with 0o Q (r) the ground state valence quark wavefunction for a quark of flavor a, with upper 
and lower components u a (r) and ia ■ rv a (r)/r, k± = u>i ± Mx, with Ui the ground state 
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struck quark energy eigenvalue, EB(r) = < EB, Rcm = r\EB, Rcm = > is the overlap 
function for two "empty bags" separated by a distance r, which accounts for the dynamics 
of the confining degrees of freedom. Finally, the subscripts si and s2 denote the flavor of 
the two spectator valence quarks that make up the nucleon. 

All calculations described in this paper are carried out using the Los Alamos Model 
Potential (LAMP)[|1^], in which valence quarks are confined by a linear potential of the form 

V{r) = /3k(r - r ) (4) 

with parameters k=0.9 GeV and r =0.57 fm chosen to reproduce the average nucleon-delta 
mass, and (3 is a Dirac gamma matrix. We further assume that the function EB(r) is a 
constant jl3H, and unless otherwise noted, quarks are taken to be massless. 



3 Charge Symmetry Breaking 

As we have already noted, the sources of charge symmetry breaking are light quark mass dif- 
ferences and the electro-magnetic interaction. CSB effects may be manifested either directly, 
as a result of explicit mass or interaction terms in the quark model hamiltonian, or indirectly 
as a result of mixing between the symmetry violating terms and the strong interaction. In 
this section, we calculate the direct terms, deferring discussion of the numerically smaller 
mixed interactions for later. 

At the partonic level, charge symmetry predicts that the u quark distribution in the 
proton is equal to the d quark distribution in the neutron, with a corresponding prediction 
for the d distribution. A measure of the extent to which the symmetry is broken is given by 
the ratios 



I'( <Sb( X ) — r /• 



2 KW-<(i)] 

Uy{x) + dy{x)\ 



r>min ( \ _ 2[(jy(x) Uy{x)\ . . 

kcsb[x) - [<(*)+<(*)], (5) 



where u{d)y n \x) denotes the up(down) valence quark distribution in the proton (neutron). 
3.1 CSB in Quark Wavefunctions 

The simplest mechanism for altering the shape of the valence distributions of the nucleon 
in Eq. 2 is to change the model wavefunctions for the struck and spectator quarks. Since 
these wavefunctions are generated by solving a Dirac equation in the mean confining field 
of the other quarks, direct CSB effects may be included without additional assumption by 
including appropriate terms in the Dirac Hamiltonian. 

We have recalculated the valence quark wavefunctions assuming quark masses m u = 4 



MeV and md= 8 MeV[14|. In Figs. 1 and 2, the CSB ratios produced by these wavefunctions 



are indicated by the dashed curves. For minority quarks, the ratio starts small and rises to 
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+1.5% at large x, while the majority quark ratio rises more slowly, but is still positive. The 
fact that these ratios are of the same sign and comparable in magnitude may be understood, 
since, when the struck quark is a majority quark, the two spectators are in a charge symmetric 
state. The only change is due to the change in the valence quark wavefunction. When a 
minority quark is struck, the change in the wavefunction of the struck quark is compensated 
for by the change in the wavefunctions of the two spectators, which alter the momentum 
projection factor Gi(k) in Eq. 1, and are of opposite sign and roughly twice as large as the 
effect produced by the struck quark wavefunction. 

In addition, we have calculated the effect of the Coulomb interaction between the charged 
quarks using a mean field approximation for the electric potential between a quark of charge 
q and the other valence quarks, given by 

V Co ul(r) = aq(Q N - q) [ dV ^'^f ) , (6) 

J |r — r | 

where Qn is the charge of the nucleon being studied. The CSB contribution generated by 
the Coulomb effect on the quark wavefunctions is shown by the dot-dashed curves in Figs. 
1 and 2. Again, the CSB effect produced by perturbing the wavefunctions is significantly 
smaller for the minority quark distributions. In this instance, the change in the wavefunction 
due to the Coulomb force is, to first order in a, the same for either minority quark. Even for 
the majority quarks, the Coulomb correction is much smaller than 1%. Magnetic corrections 
have not been calculated explicitly, but since they are similar in structure to the color 
magnetic corrections responsible for breaking SU(4) symmetry in the quark distributions 0, 
yet are suppressed relative to those corrections by both the smallness of the electromagnetic 
coupling constant and by color SU(3) factors, their effect should also be very small. 

3.2 CSB in the Dirac Eigenvalue 

Along with the wavefunctions, solution of Dirac equation provides an energy eigenvalue which 
determines the time dependence of the lowest mode of the confined quark field. This energy 
is just that required to break the bonds which tie the struck quark to the spectators, and as 
such is sensitive to the nature of the interactions that bind the quarks together. Since the 
dependence of the quark distribution functions on this eigenvalue has the same functional 
form regardless of whether the changes in the eigenvalue are produced by quark mass effects 
or Coulomb interactions, we have simply taken the changes in the eigenvalues produced by 
the Hamiltonian changes already described and added them together. For the proton the 
quark eigenvalues were shifted from the value for massless quarks (361.8 MeV) to uo u = 364.2 
MeV and Ud = 365.0 MeV. The corresponding values for the neutron are u u = 363.0 MeV 
and uid = 365.6 MeV. The CSB ratios obtained using these eigenvalues, and massless quark 
wavefunctions, are shown by the solid lines in Figs. 1 and 2. At large x, these shifts produce 
CSB effects on the order of 2-3%, which, may be understood from Eq. 1 as the fractional 
shift in the quark eigenvalue enhanced by the slope of the unperturbed quark distribution. 
(This is the closest analog to the diquark mass shift of Ref. 6, and correspondingly produces 
the largest effect.) 
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3.3 Proton-Neutron Mass Difference 

The results of changing the nucleon mass parameter appearing in Eq. 1 to reflect the 
difference in mass of the proton and neutron is shown by the dotted curves in Figs. 1 and 2. 
Like the quark eigenvalue, the effect of changing the nucleon mass parameter is enhanced, 
at large XBj, by the slope of the unperturbed quark distribution function. 

This effect, however, will not be present if the data is analyzed in the conventional 
manner because the nucleon mass appears not only as an explicit parameter in Eq. 1, 
but is also implicit in the definition of xsj (In deep inelastic leptoproduction, for example, 
XBj = Q 2 /2M nuc q in the target rest frame.). In fact, the combination M nuc XBj, which is 
all that appears in Eq. 2 for xq(x), is completely independent of M nuc \ This is in accord 
with what one would expect in a light cone formalism, where P + is a kinematical variable, 
and therefore immune to the dynamical effects which violate charge symmetry. (While it 
is certainly possible to alter the usual analysis by rescaling xsj once the neutron structure 
function is extracted from the raw data, it not clear what would be learned by comparing 
the probability of finding quarks in the proton at one momentum to that in the neutron at 
a slightly different momentum.) Hence, we assert that there is no CSB effect in the parton 
distributions due to the neutron-proton mass difference f p~5f| . 



4 CSB in the Confining Interaction 

The existence of CSB interactions induced by mixing the electromagnetic and strong cou- 
plings was first pointed out in Ref. 16, where it was argued that the quark-gluon vertex picks 
up a charge asymmetric contribution from photon loops|16|]. Since the confining potential 



must, at some level, be composed of multiple gluon exchanges fll?!, the existence of this cou- 
pling implies that the confining potential will not be charge symmetric. Unfortunately, since 
the confining potential is a Lorentz scalar, the CSB contribution to the energy of light quark 
hadrons is indistinguishable from the effect of quark mass differences. 

In this section, we examine the possibility of distinguishing CSB in the confining potential 
from quark mass effects by looking at the relative contribution of each to CSB in the valence 
quark distribution. Since there is, at present, no means to calculate the confining potential, 
nor its correction due to charge symmetry violation, we proceed to model the effect by 
altering the string tension parameter used in the LAMP model potential. Furthermore, 
since the relative normalization of the two effects is unknown, we proceed by arbitrarily 
normalizing the shift in the string tension such that it produces the same first order shift 
in the quark eigenvalue as the corresponding quark mass. The wavefunctions that result 
from this change are then used in Eq. 1 to produce the the charge symmetry breaking ratios 
shown in Fig. 3. The solid and short-dashed curves indicate the CSB ratios for majority 
and minority quarks produced by a 4 MeV quark mass difference, while the long-dashed 
and dot-dashed curves are the same ratios produced by the change in the string tension. 
Significantly, the CSB produced by the change in the string tension has the opposite sign to 
that produced by quark mass differences, opening the possibility that the two effects may 
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be distinguished from one another by precise measurement of CSB in the valence quark 
distributions. 



5 Evolution 

Having separated out the CSB effects at a low momentum scale, we must now evolve to high 
Q 2 so that a comparison with data is possible. Two issues arise: First, does the perturbative 
evolution erase the CSB effects we have calculated, and secondly, how large are the CSB 
effects in the evolution itself? 

In Figs. 4 and 5, we show the results obtained by evolving both the CSB pieces of the 
valence quark distributions and the symmetric distributions from a low quark model scale, 
taken to be 0.5 GeV 2 , to 10 GeV 2 ||. For both majority and minority quarks the evolution 
has only a small effect on the CSB ratios, shifting the curves to low x and slightly increasing 
the magnitude of the ratio. 

Also shown in Figs. 4 and 5 is the contribution to the CSB ratio provided by charge 
symmetry breaking effects generated when a quark splits into a quark and a photon [|I~8||. 
To leading order in a, the structure of the QED and QCD contributions to the evolution 
differ only by constant factors, and the effect of including the photon diagram is to slightly 
speed up the rate at which the valence distributions evolve. Since the coupling of the photon 
is charge asymmetric, this means that the u distribution evolves slightly faster than the d 
distribution. As shown in the figures, this effect is larger than the effect of Coulomb repulsion 
in the quark wavefunctions, but is nonetheless still quite small. 



6 Comparison and Experimental Prospects 

In Ref. 6, the effect of charge symmetry violation on the valence quark distributions was 
calculated via the introduction of an intermediate state diquark. The resulting two-body 
kinematics produces an additional enhancement of the CSB effect produced by small changes 
in the diquark mass, resulting in a 5-10% CSB ratio for minority quarks in the range 0.5 < 
x < 0.7, roughly twice the size of our result. For majority quarks, the two calculations are 
comparable. In reference 5, using a model independent approach, Sather obtains slightly 
smaller ratios for the majority quarks than in this work, and similar results for the minority 
quarks. Generally, each of these calculations predict larger CSB effects in the minority quark 
distributions than in the majority distributions. 

Experimentally, the situation is less clear. As pointed out in Ref. 18, measurements of 
CSB generally yield not the ratios for the majority /minority quarks separately, but rather 
the ratio for the entire valence distribution, given by 

2[d v {x) - u v {x) - u v {x) + d v {x)} 

^ [U P V (X) + d P y(x) + Uy(x) + dy(x)} ' ^ ) 
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Since the minority quark distribution is suppressed at large XBj by SU(6) symmetry breaking 
effects [ p0|1 , the ratio is less sensitive to the (fractionally) larger charge symmetry breaking 



in the minority distribution, which is the major difference between the models. While it is 
possible to isolate the minority distribution by comparing tt + — p and n~ — D Drell-Yan cross 
sections [|19|, the systematic errors associated with beam normalizations and the deuteron 



EMC effect will make it difficult to unambiguously separate the small effects predicted in 
this work. 

Regardless of the result, the information provided by these experiments will provide new 
insights into the soft dynamics of quarks in the nucleon, and possibly a means to distinguish 
experimentally between CSB generated by quark mass differences and the mixing of the 
strong and electro-magnetic interactions. 
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Figure Captions 



Figure 1 - Charge Symmetry Ratio for Majority (valence) quarks in the nucleon at the 
quark model scale. 

Figure 2 - Charge Symmetry Ratio for minority (valence) quarks in the nucleon at the 
quark model scale. 

Figure 3 - Comparison of Charge Symmetry breaking produced by quark mass differ- 
ences with that induced by mixing of the electro-magnetic and confining interactions. 

Figure 4 - Charge Symmetry Ratio for Majority (valence) quarks in the nucleon at a 
scale of 10 GeV 2 . The sum of all contributions is shown by the solid line, while the 
other lines are described in the text. 

Figure 5 - Charge Symmetry Ratio for Minority (valence) quarks in the nucleon at a 
scale of 10 GeV 2 . The sum of all contributions is shown by the solid line, while the 
other lines are described in the text. 
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